Emergent Geometric Hamiltonian and Insulator-Superfluid Phase Transitions 
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I argue that certain bosonic insulator-superfluid phase transitions as an interaction constant 
varies are driven by emergent geometric properties of insulating states. The renormalized chemical 
potential and population of disordered bosons at different levels define the geometric aspect of an 
effective low energy Hamiltonian which 1 employ to study various resonating states and quantum 
phase transitions. In a mean field approximation, I also demonstrate that the quantum phase 
transitions are in the universality class of a percolation problem. 
PACS number: 03.75 Kk, 03.75.Lm, 03.70.+K 
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Localization or derealization of interacting bosons has 
been a puzzling theoretical and experimental issue for 
quite long timei*2iMiSiSiL2iS. Anderson et al. first in- 
vestigated a possibility of bose condensation into an ex- 
tended eigenstate in disordered systems and pointed out 
certain instabilities of such condensates^. Giamarchi 
and Schulz employed the replica-trick-based renormal- 
ization equations to study the quantum phase transi- 
tions between superfluid states and insulating states in 
one dimension^. Fisher et al. on the other hand in- 
vestigated the localization of interacting bosons in the 
context of a Bose-Hubbard (lattice) model&. In high- 
dimensions, this different approach appears to be very 
fruitful resulting in various new observations; for in- 
stance a new distinct Bose-glass phase has been em- 
phasized and its divergent superfluid susceptibility was 
pointed out. In the same context, numerical results 
further suggested additional Anderson-type insulating 
states in weak couplings^. More recently, Id supcrfluid- 
Mott-glass phase transitions in a particle-hole symmetric 
lattice model have been studied 9 . Finally, a system of 
cold atoms in a random laser speckle potential has also 
been established and studied^. 

However, little attention has been so far paid to the 
microscopies of bosonic insulators in the dilute limit. 
Technically, the complexity of this issue is largely related 
to the failure of the replica-technique-based perturbative 
approach applied in this limit (in both ID and high di- 
mensions). This aspect of the problem seems likely to 
originate from a very high degree of inhomogeneity of 
bosonic insulators which is hard for a standard pertur- 
bative replica-trick-based approach (either mean field or 
one-loop expansion) to capture. When all low energy 
one-particle states are localized, the instability of a con- 
ventional ideal boson condensate (with a given number 
density in an infinity system) towards a local condensate 
with infinite number density appears to be evident. It is 
the finite repulsive interaction which prevents such col- 
lapsing from taking place. This distinctly differs from the 
problem of fermion localization: various fixed points in 
the one-parameter scaling theory for noninteracting sys- 
tem seem to be stable in weakly interacting systems^. 

In this Letter, I would like to study the renormalized 
chemical potential and population of bosons at differ- 



ent localized states to illustrate the microscopic struc- 
ture of bosonic insulators. Secondly, I will argue that 
interactions between bosons result in an emergent ge- 
ometry of bosonic insulators. This observation suggests 
that certain superfluid-insulator quantum phase transi- 
tions should be predominantly driven by geometrically 
percolating localized states. I also apply the mean field 
solutions to probe the low energy Hilbert subspace rel- 
evant to quantum critical points. At last I propose to 
study critical behaviors of these quantum phase transi- 
tions using an effective Hamiltonian defined in percolat- 
ing clusters. 

Consider bosons in a d-dimension space with simple 
short range random impurity potentials. All relevant low 
energy one-particle states are already localized (and are 
far away from the mobility edge which separates extended 
one-particle states and localized states). Furthermore, 
bosons interact via a weak two-body delta-like interac- 
tion UoS(r — r'). The set of wave functions I am going to 
exam can be written in the following form 
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bt is the creation operator of a boson at a localized one- 
particle state k of energy and 4>k is the corresponding 
phase. Different states correspond to different choices of 
g({rik})- For an insulator in the weakly interacting limit 
(i.e. the dimensionless interaction constant T is small), 
g({rik}) is nonzero only for a given set of integer values 
of rik = Int[n(ek)] and no resonance is involved. n(e) as 
a function of e is calculated below and depends on both 
interactions and disorder strength. Here Int[I + rj] takes 
integer value I if 1/2 > r] > and I + 1 if 1 > 77 > 1/2; 
at 77 = 1/2 it is equal to either I or I + 1/2. For a 
superfluid phase, g is a smooth function within a wide 
range of distribution {n^}; an extreme case is the equal- 
amplitude state with g({rik}) being constant. 

Consider an insulator in the weakly interacting limit 
where the spatial overlap between different occupied 
states is negligible; g({n k }) = ]J k 5 nk jnt[n(t k )}- The total 
(mean field) energy for such a state which is fully charac- 
terized by {n(efc)} and is independent of {4>k} (see Eq.JQ 
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and discussions following it) can be shown to be 



E ({n(e k )},{M) = J de{e-p)v{e)n{e) 
+U Jdede'n(e)n(e')C(e,e'), (2) 
C(e,e') = J d d vv{e,v)v{e' ,v). 

p is the chemical potential to be determined self- 
consistently. z/(e) and i/(e, r) are the density of states and 
local density of states respectively, ^(e,r) = J2k^( e ~ 
ek)tpl(r)ipk{r), v{t) = fd d rv(e,r). k is a quantum 
number labeling one-particle eigenstates in the disor- 
dered system. Finally, the distribution is subject to a 
constraint that the average number density is po, i.e., 
J den(e)v(e)— Vopo- (Vo is the Volume.) Minimizing 
E ({n(ek)}) with respect to {n(efc)} suggests the follow- 
ing integral equation for the ground state population dis- 
tribution (e < p) 



2U / de'n(e')C(e, e') = v{e){p - e). 



(3) 



I practically seek for a solution when interactions are 
weak and the chemical potential p is much smaller than 
A = (^(m)Cl) 1 ! the typical single particle level spacing 
at chemical potential in a volume defined by the localiza- 
tion length £l (And for simplicity, we neglect the energy 
dependence of £l for all states under consideration are 
far away from the mobility edge.). In the leading order 
of /Li/ Ao, I obtain the solution to Eq. © (after an average 
over random potentials) in d-dimensioni^ 



n(e) 



N L p 



2r 
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when e < p; otherwise, n(e > p) = 0. Two dimensionless 
parameters entering the solution are 



(5) 



which characterize the strength of disorder and interac- 
tions respectively, cf — p 2 J d jlm is the Fermi energy 
for free fermions with the identical mass (m) and density 
(po); S is a constant of order of unity (I have set K = 1). I 
also assume that Nl is much larger than one. It is worth 
remarking that p, the renormalized chemical potential of 
disordered bosons is proportional to pg^ d+2 \ 

The microscopic properties of the insulating states are 
fully characterized by this distribution function. For in- 
stance, one can evaluate £xj the typical distance between 
two states occupied by weakly interacting bosons, the 
typical number density of bosons in each occupied state 
Pl and the ratio between Aq and p. I list the results here 



Ao 
A* 



Pa 



N, 



N Ljd/(d+2) 
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which have been expressed in terms of Nl and T. Follow- 
ing Eq.© that in a dilute limit where £t is much larger 

than £l, or r <C N l 2 ^ d (<^ 1), the local level spacing Ao 
is much larger than //. 

If the interaction is indeed very weak so that £t ^ £l, 
then typical finite-size condensates at various localized 
states with the density pl are spatially disconnected. 
For a typical localized state k, the number of particles 
in the ground state is unique except when n(ek) is a half- 
integer. The energy gap in local particle- like or hole-like 
excitation spectra can be shown to be proportional to 
{1 ± 2Int[n(ek)] T 2n(efe)} Ypf/Nl and distribute in a 
range between to aoTpp/NL (ao is a quantity of order 
of unity) . One therefore obtains the typical on-site tem- 
poral correlation function G T (r, 0) =< 6 + (t, r)6(0, r) >, 



G T (T,0Hexp(--),T c ~-^ 



(7) 



which reflects the well-known J7(l)-symmetry restoring 
occurring independently at different localized states, (the 
superscript T indicates a typical value). For the non- 
interacting case, r c goes to infinity as it should for ideal 
bosons. 

When n(efe) is precisely a half integer, however low- 
est energy states are two-fold degenerate. These states 
play an important role in the calculation of the averaged 
Green's function and can yield to an ensemble averaged 
power-law temporal correlation and divergent suscepti- 
bility. 

Eq.0} indicates that locally, many-body collective 
states with Int[n(ek)] Int[n(ek)] ± 1-, Int[n(ek)] ± 2-, ... 
bosons form a low energy Hilbert subspace. To establish 
a phase coherence and superfluidity, a local condensate 
k with nk — Int[n(ek)] bosons has to resonate with low 
lying states of nk ± 1-, nk ± 2-, ... bosons. For stud- 
ies of the physics close to insulator-Superfluid quantum 
phase transitions, one has to therefore take into account 
all the low energy states represented by = Int[n(ek)], 
Int[n(ek)] ± 1, Int[n(ek)] ± 2 etc. A convenient basis to 
discuss resonance between these states is the following 
coherent state representation 



\{M >~ 11zJ-9oK - Int[n(e k )}) k ' \vac >, 
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g (n) is unity for n max > \n\ > 0; n max is much larger 
than one but smaller than n(ek) so that all the low energy 
states to be involved in resonance are present in this 
representation. 
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The energy of these states, in addition to that in 
Eq. J2Jl acquires contributions from the the exchange 
(random) coupling between two resonating condensates. 
This consideration and Eqs.©, © indicate that close 
to insulator-superfluid transitions, the effective Hamil- 
tonian in the truncated low energy subspace spanned by 
states in Eq.© should bei^ 



Heff = Y,E c (k)(fi k -n(e k )) 2 
allk 

Jki cos(0 fc - (pi)-, 

C <fcZ>£C 



(9) 



Here h k — b\pk is the number operator of bosons at state 
k; and [hk,4>i] — —iSk.k'- The magnitude of n(efc) is 
given in Eq.Q. The amplitude of Jki is given by J = 
Pl^lM; B k ~ if J dr|^(r)V'fc(r)| 2 , and the overlap inte- 
gral between two localized states is Aki ~ / ipl(r)ij}i(r) 
p{y)/pl<Iy] p{r) is the local boson density. Finally, in 
deriving Eq. © I have neglected a) exchange interactions 
between pairs of states of distance much longer than £l 
because Aki , the overlap integrals of those pairs are expo- 
nentially small; b) exchange couplings involving four con- 
tacting localized condensates. Only exchange couplings 
between two contacting localized states separated by dis- 
tance smaller than £l (labeled as < kl > in the second 
sum) are kept. Configurations of contacting localized 
states define clusters {C} of various sizes and represent 
resonating localized condensates. Note that for a given 
r, properties of an insulating state are determined by the 
Hamiltonian in Eq.© defined in percolating clusters {C} 
of resonating states; the distribution of clusters {C} itself 
is a function of T as an emergent property of insulating 
states. 

The ratio between the above two energies, l/r c and J 



1 = „( r \2d/(d+2) 



(10) 



which is much less than unity when Y <C Nl- (a and b 
introduced below are constants of order of unity.) Close 
to the percolation threshold discussed below, 1/Jt c is 
proportional to 1/N 2 and turns out to be much smaller 
than unity as Nl in our case is much less than one. 

In a mean field approximation where 1 / Jr c is set to 
be zero, phases of two localized condensates are locked 
once they are in contact. Gj(r, 0), the typical temporal 
Green's function of bosons belonging to a given cluster C 
(with random chemical potentials of course), is a function 
of the size s of the cluster C only. For a given cluster, this 
Green's function is still given by Eq.©; however, t c is 
replaced with st c implying the slow down of the quantum 
symmetry restoring in large clusters. 



Solutions to Eq.© can be obtained in this mean field 
approximation (d > 2 for the rest of discussions). When 
all clusters {Coo} present are infinite, the following wave 
function of a resonating state with various clusters rep- 
resents one of such solutions^, 



t p i<t>k\rik 



(Pie 



Coo fceCoo n k 



IIE Yl — [ —\v<><- 



(11) 



where the second term corresponds to the wave function 
of localized bosons that do not belong to any clusters 
Coo. The set of variables {(f>k,k € Coo} minimize the 
exchange couplings in Eq.©. Again go(nk — Int[n{ek)]) 
is a smooth function in the vicinity of nk = Int[n(ek)] 
which stands for resonance between localized condensates 
with different occupation numbers. 

Close to the percolation threshold, one obtains the fol- 
lowing formula of the Green's function 



G T (r,0) = ]TP( S )exp(- — ); 



(12) 



P(s) is the probability of finding a localized state be- 
longing to a cluster of size s and exhibits interes ting 
scaling behaviors (for general discussions, see RefUJ, 
Il3^ . Let us define that the concentration of occupied 
states as p(T,N L ) = b{£, L /£, T ) d = bN L {T / N L ) d ^ d+2 \ 
The one-parameter scaling ansatz suggests that P(s) ~ 
Po(s/£(p)), Thus, one arrives at the following scaling 
form G t (t, 0) ~ g t {r / V c £(p)); gt{x) approaches zero 
when x goes to infinity. Here £ (p) is the percolation cor- 
relation length expected to diverge as \p — p c \~ u ', close 
to the critical interacting constant , one generally has 

£(r,Jv L )~(r£-r)-\ 

Ea. l|12|l and discussions above therefore indicate that 
t/(l)-symmetry breaking should take place if and only if 
the percolating threshold is reached so that £(p) is in- 
finite. The appearance of infinitely percolating clusters 
(with a finite density) turns out to be necessary and suf- 
ficient (in a mean field approximation) for the develop- 
ment of long range correlations. An estimate based on 
Eq.© yields the critical point for geometric quantum 
phase transitions 



Y L = (Plyd/{d+2) N -2/d 



(13) 



The superscript L stands for a lower critical point for the 
superfluid phase. As d approaches infinity, the percola- 
tion threshold p c scales as 1 / d and clusters of infinite size 
first appear in the dilute limit where £t 3> £l- The per- 
colation can be studied in a mean field approximation 15 . 

Furthermore, the superfluid density above the perco- 
lation threshold p c can also be related to clusters ex- 
tending to infinity. One can easily establish that , 
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the upper bound of the superfluid density p s should be 
proportional to n^fT), the probability density for a lo- 
cal state to be in an infinite cluster appearing above T^. 
P U S = Po«-oo(r) ~ (T ~T^)^. In the mean field approxi- 
mation, the critical exponents are v = 1/2, j3 = 1. 

The above mean field approximation appears to be 
oversimplified from the point of view of Shklovskii- 
De Gennes backbone (SDeGb) model for percolating 
cluster a 16 ! 17 . In fact, the SDeGb model indicates that 
percolating clusters consist of dense regions connected by 
one-dimensional strands even when d > 1. Early results 
further imply that it should be unlikely for a ID infi- 
nite quantum system to be ordered^. So at any small 
but finite value of 1/Jr c , we expect, after taking into 
account quantum effects, that the onset superfluidity ap- 
pears above the percolation threshold. Particularly, the 
true quantum transition takes place when bosons along 
ID strands (with finite length) in the SDeGb model be- 
come correlated in phases. Thus, the lower critical value 
r^g receives a quantum correction, 



1 cQ 



6d(l/Jr c ) 



l/ 7 



(14) 



Here £,id(&) is the correlation length of bosons along 
ID strands in a Shklovskii-de Gennes cluster and is ex- 
pected to be exponentially long compared to £l when 
a approachs zero. 7 is the critical exponent of percola- 
tion correlation length above the threshold. As indicated 
in Eq. l|14|l . the quantum correction to the lower critical 
point r^g — is negligible when 1 / Jt c is much smaller 
than unity. Though in this case the critical phenomenon 
might differ from the percolative ones, the size of the 



critical regime where the deviation from the mean field 
critical behavior becomes evident could be exponentially 
small as Ea. l|14|l which sets this size implies. 

At last let us consider what happens when interactions 
are further increased and T ^> N^. Bosons at different 
localized states then prefer to have independent dynam- 
ics; that is, Gc(t, 0;0,0) ~ exp(— t/t c ) and therefore 
G T (r,0) discussed in Ea. lpH) remains short ranged even 
when there are infinite clusters^. The above argument 
suggests that the upper critical point for the superfluid 
phase should be set by the following equation: = cNl 
and c is a number of order of unity. The superfluid phase 
appears when > r > . 

When Nl is much less than one as in the presence 
of strong disorder, <g; and the regime where the 
superfluid phase exists shrinks to zero. Indeed, bosons in 
each localized state have independent dynamics, either 
because of the absence of connected clusters (T <C T^) 
or when above the percolation threshold, because of the 
rapid quantum symmetry restoring (T > 3> ). No 
insulator-superfluid phase transitions are likely to occur 
in this case and there is one single insulating phased. 

To summarize, I show that bosonic insulators of in- 
teracting bosons have fascinating geometric properties 
which are characterized by the parameters in Eq.©. 
Certain quantum phase transitions from insulators to 
superfluid phases as the interaction constant varies are 
strongly influenced by the emergent geometries of insu- 
lating states. This work is in part supported by NSERC, 
Canada and a grant from UBC. I would like to thank 
P. W. Anderson, B. Bergersen, G. Refael, S. Sachdev, 
B.Spivak and P. B. Wiegmann for helpful discussions. 



1 P. W. Anderson, Mater. Res. Bull. 5, 549 (1970); J. A. 
Hertz, L. Fleishman, and P. W. Anderson, Phys. Rev. Lett. 
43, 942 (1979). 

2 J. D. Reppy, Physica 126B, 335 (1984); J. Low Temp. 
Phys. 87, 205 (1992). 

3 M. Ma, B. I. Halperin, and P. A. Lee, Phys. Rev. B34, 
3136 (1986). 

4 T. Giamarchi and H. J. Schulz, Phys. Rev. B37, 325 
(1988). 

5 D. Finotello, K. A. Gillis, A. Wong, and M. H. W. Chan, 
Phys. Rev. Lett. 61, 1954 (1988). 

6 M. P. Fisher, P.B. Weichman, G. Grinstein and D. Fisher, 
Phys. Rev. B40, 546 (1989). 

7 R. T. Scalettar, G. G. Batrouni, and G. Zimanyi, Phys. 
Rev. Lett. 66, 3144 (1991). 

8 W. Krauth, N. Trivedi, and D. Ceperley, Phys. Rev. 
Lett. 67, 2307 (1991). 

9 E. Altman, Y. Kafri, A. Polkovnikov, and G. Refael, Phys. 
Rev. Lett. 93, 150402 (2004). 

J. E. Lye, L. Fallani, M. Modugno, D. Wiersma, C. Fort 
and M. Inguscio, cond-mat/0412167 See also R. Roth and 
K. Burnett, Phys. Rev. A 68, 023604 (2003); B. Damski, 
J. Zakrzewski, L. Santos, P. Zoller, and M. Lewenstein, 



Phys. Rev. Lett. 91, 080403 (2003). 

1 P. W. Anderson, Phys. Rev. 109, 1492 (1958); D. J. Thou- 
less, Phys. Rev. Lett. 39, 1167 (1977); E. Abrahams, P. 
W. Anderson, D. C. Licciardello, and T. V. Ramakrishnan, 
Phys. Rev. Lett. 42, 673 (1979); L. P. Gorkov, A. I. Larkin 
and D. E. Khmelnitsky, JETP Lett. 30, 228 (1979); See 
also P. A. Lee and T. V. Ramakrishnan, Rev. Mod. Phys. 
57, 287 (1985). 

2 M. Plischke and B. Bergersen, Equilibrium Statistical Me- 
chanics (World Scientific, 1994). 

3 P. M. Chaikin and T. C. Lubensky, Principles of condensed 
matter physics, (Cambridge Press, 2000). 

4 F. Zhou, unpublished (2005). 

5 The effect of level repulsion on percolation can be taken 
into account by redefining the concentration p. 

6 A. S. Skal and B. I. Shklovskii, Sov. Phys. Semicond. 8, 
1029 (1975). 

7 P. G. De Gennes, J. Phys. Paris. Lett. 37, LI (1976). 

8 J. W. Hartman and P. B. Weichman, Phys. Rev. Lett. 74, 
4584 (1995). 

9 The impurity potentials considered in this article are un- 
bounded; for bound ed potentials as those introduced and 
studied in Ref l6l7l8l there can be more than one insulating 



■5 



phases. 

For cold atoms, the dimensionless coupling constant T can 
be experimentally varied either by varying the two-body 



scattering length close to the Feshbach resonance or simply 
varying the density po- 



